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BEAM  CURRENT  LIMITATIONS  DUE  TO  INSTABILITIES  IN 
MODIFIED  AND  CONVENTIONAL  BETATRONS 

I.  Introduction 

The  addition  of  a  toroidal  magnetic  field  to  a  conventional  betatron  has 

been  shown  theoretically^-^  to  increase  the  equilibrium  current  which  may  be 

2 

confined  by  a  factor  1/2  (B  /B  )  ,  for  large  values  of  B  /B  .  In  conventional 

0  Z  0  z 

accelerators  without  solenoidal  focusing,  however,  beam  stability 

6—8 

considerations  place  the  actual  limit  on  beam  current.  Therefore  it 
becomes  important  to  analyze  the  stability  conditions  and  associated  limiting 
currents  for  a  given  beam  equilibrium  in  the  presence  of  a  toroidal  field.  In 
this  paper  we  present  such  an  analysis  for  both  longitudinal  and  transverse 
modes. 

A  device  in  which  a  toroidal  magnetic  field  is  superimposed  on  the  usual 
weak  focusing  betatron  field  has  come  to  be  called  a  "Modified  Betatron".  See 
Fig.  1.  A  stability  analysis  of  this  accelerator  necessarily  must  include  the 
strong  self  and  induced  (wall  image)  fields  of  the  electron  beam.  It  is 
primarily  the  inclusion  of  these  fields  which  distinguishes  this  work  from  the 

9 

stability  analysis  performed  for  the  so-called  plasma  betatron  in  which  self 
fields  are  much  less  important.  These  self  field  effects,  however,  will  be 
seen  to  have  a  dramatic  effect  on  the  current  versus  energy  spread  scaling; 
namely  we  predict  the  existence  of  more  than  one  stable  value  of  current  for  a 
given  beam  energy  spread.  This  somewhat  surprising  result  will  be  discussed 
later.  The  central  result  of  this  work,  however,  is  that  significantly  more 
current  may  be  carried  by  a  beam  in  a  modified  betatron  configuration  than  in 
a  conventional  betatron  of  the  same  dimensions,  assuming  equal  beam  sizes  and 
energy  spreads.  In  reaching  this  conclusion  we  have  included  not  only  the 
effects  of  self  fields,  but  also  the  (stabilizing)  effects  of  betatron 
oscillations  and  the  (destabilizing)  effects  of  short  wavelength  enhancements 
to  the  longitudinal  and  transverse  impedances  due  to  chamber  resonances. 
Manuscript  approved  July  12, 1983. 


Below  we  discuss  the  dispersion  relation  for  arbitrary  toroidal  fields  and 
currents,  describe  our  model  for  the  impedances,  and  present  analytical  and 
numerical  results  from  the  dispersion  relation. 

II.  Discussion 

A  dispersion  relationship  for  both  the  longitudinal  and  transverse  modes 
in  a  modified  betatron  accelerator  configuration  has  been  derived. ^  Included 
in  the  derivation  are:  beam  self  field  effects,  induced  field  effects  arising 
from  wall  image  charges  and  currents  as  well  as  finite  chamber  wall 
conductivity  effects.  Toroidal  corrections  to  the  equilibrium  beam  self 
fields  and  chamber  wall  image  fields  have  been  neglected.  The  longitudinal 
and  transverse  impedances,  which  characterize  the  beam  environment,  will  here 
be  incorporated  in  a  phenomenological  way  in  the  short  wavelength  limit.  The 
dispersion  relation,  therefore,  treats  disturbances  of  all  wavelengths. 
Including  wavelengths  much  longer  or  shorter  than  the  chamber  minor  radius. 

The  short  wavelength  model  for  the  impedances  contains  effects  associated 
with  propagating  chamber  modes.  These  effects  can  significantly  affect  the 
instability  growth  rates.  With  the  inclusion  of  the  short  wavelength 
contributions  to  the  impedances  a  realistic  estimate  for  the  current 
limitations  due  to  the  various  instabilities  can  be  obtained  for  the  modified 
betatron  configuration  and  compared  to  those  of  a  conventional  betatron.  To 
perform  a  meaningful  comparison  we  will  choose  identical  parameters  for  the 
two  types  of  betatrons,  i.e.,  same  geometry,  injection  energy,  field  index, 
etc.  The  only  difference  of  course  will  be  that  the  modified  betatron 
configuration  will  include  a  toroidal  magnetic  field. 

The  dispersion  relation^-®  for  the  longitudinal  and  transverse  modes  of  a 


cold  beam  may  be  written 


where  Aw  *  ( a>  -  Zw  )/w  ,  w  is  the  complex  mode  frequency,  Z  *  1,2, 3,... is  the 
Zee 

longitudinal  (toroidal)  harmonic  mode  number,  u>  *  ft  /y  is  the  electron 

c  z 

rotation  frequency,  ft  *  |e|B  /m  c  is  the  non-relativistic  cyclotron 

z  z  o 

2  2  -1/2 

frequency,  y  ■  (1  -  v  /c  )  ,  vq  is  the  longitudinal  (toroidal)  beam 

^  2  2 

velocity,  w  *  -  2iZ  (v/y)(Z  ^/zZq),  u  is  Budker's  parameter 

3 

(v  ■  ITA1/17  x  10  ),  Z  ■  Z  (m)  is  the  total  effective  longitudinal  impedance, 

i!  n 

ZQ  -  Wc  (ZQ  -  377  a  in  MKS  units),  Wj2-  1/2 
3  2 

-  2(v/y  )(r  /r,)  -  2i(v/y)(r  Z,/Z  ),  r  is  the  major  electron  beam  radius, 

O  D  OXOO 

r^  is  the  minor  electron  beam  radius,  Z^  »  Z^(u>)  *s  t'16  tota^  effective 

transverse  impedance,  b  ■  B  /B  and  B  is  the  toroidal  magnetic  field.  In  (1) 

9  z  9 

finite  amplitude  betatron  oscillations  were  neglected,  the  external  field 
index  was  taken  to  be  1/2  and  the  electron  beam  was  assumed  to  be  mono- 
energetic,  highly  relativistic  (vq-  c)  and  circular  in  cross  section. 

III.  Approximate  Representation  of  Impedances 

In  our  model  the  longitudinal  impedance  ^(w)  is  taken  to  consist  of 
three  terms 


Z  +  Z,  +  Z  • 
t  ,s  n  ,o  n  *  r 


(2) 


The  first  term  in  (2)  is  the  long  wavelength  space  charge  shielding 
contribution  associated  with  a  smooth  infinitely  conducting  chamber  and  is 
given  by 


Z .  -  l  Z  (i/2y  2)(1  +  2  Zn  a/r. ), 


(2a) 


where  yw  is  the  relativistic  factor  corresponding  to  the  wave  phase  velocity 

2 

and  a  is  the  minor  radius  of  the  toroidal  chamber.  Due  to  the  1/y  factor 

'w 

this  impedance  term  is  typically  quite  small;  it  may  even  change  sign. 

The  part  of  the  longitudinal  impedance  due  to  the  resistive  nature  of  the 
chamber  wall,  in  the  long  wavelength  regime,  is 

Z„  -  %  Z  (1  -  i)  6/2a.  (2b) 

a ,  0  o 

“1/2 

where  5  *  c  ( 2iro I  u>l  )  is  the  skin  depth  and  a  is  the  wall  conductivity.  It 
has  been  assumed  in  (2b)  that  the  skin  depth  is  small  compared  to  the 
thickness  of  the  chamber  wall. 

Finally,  the  last  term  in  (2)  represents  the  resonance  contribution  to 

the  total  longitudinal  impedance  and  arises  from  the  fact  that  the  chamber  can 

support  propagating  waves.  To  obtain  the  exact  form  for  Z  would  require  a 

n ,  r 

rather  involved  analysis  of  the  beam-chamber  structure  and  is  beyond  the  goals 
of  the  present  paper.  We  will,  therefore,  simply  represent  this  contribution 
by  the  phenomenological  expression1- ^ 

a)/a>r  -  i  Q(1  -  (u>/(Dr)2) 

(a»/cor)2  +  Q2(l  “  (a>/mr)2)2 

where  RQ  defines  the  chamber  shunt  impedance,  u  »  2.4  c/a  is  the  cutoff 
frequency  of  the  lowest  order  chamber  mode,  and  Q  is  the  quality  factor 
associated  with  the  chamber.  The  chamber  shunt  impedance,  RQ,  can  be 
estimated  by  noting  that  near  a  resonant  frequency  u>  *  u)r  the  longitudinal 
impedance  is  roughly  equal  to  the  free  space  impedance 

(/3  *  O  *1/3'  <3) 


z  -  r  0L.  (. 
I,r  o  w 


4 


1/3 

It  follows,  therefore,  that  RQ  *  Zq  where  i  is  the  toroidal  mode  number 

associated  with  the  resonant  frequency,  i.e.,  I r  /ujc - 

The  expression  for  the  total  transverse  impedance  is  also  written 

as  the  sum  of  three  contributions 


(4) 


The  long  wavelength  space  charge  part  of  the  impedance,  for  a  smooth 
infinitely  conducting  chamber,  is 

21,S  ‘  1  roVl/rb2  '  (4a) 

The  second  term  is  the  long  wavelength  contribution  due  to  the  resistive 
nature  of  the  wall  and  is  given  by 

Z.  -  r  Z  (1  -  i)  5/a3.  (4b) 

X,  o  o  o 

Chamber  resonances  contribute  to  the  transverse  impedance  a  part  which  we  will 
represent  by  the  form^ 


Z 


l,r 


(4c) 


where  Z.  is  defined  in  (2c). 

In  the  present  work  we  will  not  be  concerned  with  resistive  wall  effects 
which  lead,  in  conventional  acclerators,  to  well  known  longitudinal  and 
transverse  instabilities^  having,  however,  comparatively  slow  growth  rates. 


Equations  (2b)  and  (4b)  are  included  here  only  for  reference  purposes. 


We  remark  that  with  the  impedances  as  defined  above  the  dispersion 
relation  is  virtually  independent  of  the  beam  minor  radius,  r.  ,  except  for  the 

D 

weak  logarithmic  dependence  in  (2a).  This  is  reasonable  if  we  think  of  the 
dynamics  involved  in  the  various  instabilities,  that  is,  if  we  recall  that  the 
beam  centroid  moves  transversely,  even  in  the  "longitudinal"  or  negative  mass 
Instability  in  which  beam  bunching  occurs  as  the  beam  centroid  moves  in  or  out 
radially.  Motion  of  the  beam  centroid  is  affected  by  the  externally  applied 
fields,  including  those  due  to  wall  image  charges  and  currents;  these  fields, 
unlike  beam  self  fields  which  are  carried  along  transversely  by  the  moving 
beam,  do  not  depend  on  the  minor  radius  of  the  beam. 

IV.  Stability  Condition  and  Limiting  Current 

To  obtain  the  limiting  current,  based  on  stability  requirements,  for  the 
modified  and  conventional  betatron  a  stability  criterion  is  needed.  If  the 
distribution  in  particle  rotation  frequencies  is  Lorentzian  in  shape  the 
criterion  for  stability  is  simply 

T  <  tlAQl,  (5) 

where  T  is  the  growth  rate  in  the  absence  of  a  frequency  spread  and  Aft  is  the 
half  width  of  the  frequency  spread  on  the  beam.  It  should  be  noted  that  the 
large  tails  associated  with  a  Lorentzian  distribution  make  the  criterion  in 
(5)  somewhat  less  stringent  than  would  be  the  case  if  a  more  realistic  choice 
of  frequency  distribution  were  used.  However,  the  use  of  a  more  realistic 
distribution  of  particle  frequencies  would  result  in  a  considerably  more 
involved  stability  criterion.  S  ' nee  we  ar-  interested  here  mainly  in  making  a 
comparison  of  the  limiting  current'  if  tua  modified  and  conventional 
betatrons,  consistent  use  of  a  Lorentzian  for  both  devices  should  serve  our 


purpose 


Both  an  intrinsic  longitudinal  energy  spread  on  the  beam  electrons  as 
well  as  finite  amplitude  betatron  oscillations  will  produce  a  spread  in 
rotation  frequencies.  Solving  the  particle  orbit  equations  correct  to  second 
order  in  the  betatron  oscillation  amplitude,  with  self-field  effects  and 
intrinsic  energy  spread  effects  included,  the  average  spread  in  the  beam's 
rotation  frequency  is  found  to  be  given  by 


I  ah  I  -  I  +  UnlB, 


(6) 


where 


I  Aftl^  -  j  ujc  I  al  ( AE/E)  (6a) 

I  AftL  -  a)  (r. /2r  )2 1  b2/2-n  -  3/2 1  ,  (6b) 

o  C  D  0  S 

-1  -2 

with  a  *  (1/2  -  n  )  -  y  •  The  two  terms  on  the  right  hand  side  of  (6)  are 

s 

due  to  an  intrinsic  energy  spread  and  to  finite  amplitude  betatron 

oscillations,  respectively.  In  (6a)  the  fractional  intrinsic  energy  spread  i 

2  ’ 

denoted  by  aE/E  where  E  ■  ymQc  is  the  beam  particle  energy,  (y  »  1)  and 

<*  2 

ng  *  2(v/yJ)(rQ/rjj)i  is  the  field  index  associated  with  the  beam's  self 
fields.  In  obtaining  (6a, b)  we  assumed  a  circular  beam  cross  section  and  an 
external  field  index  of  1/2.  Note  that  it  is  here,  in  the  relation  between 
AE/E  and  Aft,  that  self  fields  play  an  important  role. 

By  utilizing  the  beam  envelope  equation,  the  frequency  spread 
term  lAftl.  can  be  expressed  in  a  more  illuminating  form.  The  condition  for  a 

D 

I? 

matched  beam,  i.e.  non-oscillating  minor  beam  radius,  is  “ 


(7) 


-  V  1/2  +  (r0en)  /(rb S)  . 

where  e  is  the  normalized  transverse  beam  emittance  as  measured  in  the  Larmor 
n 

frame.  Substituting  (7)  into  (6b)  gives 

I  agIb  -  ^  1(£N/rb)2/2  '  I  (rb/ro)2,f2  +  v/2^*  (8) 

y 

In  (8),  the  first  term  in  brackets  is  the  familiar  longitudinal  energy  spread 

due  to  emittance,  the  second  term  is  a  toroidal  correction  to  the  first,  while 

the  last  term  is  the  energy  spread  associated  with  the  electrostatic  potential 

drop  across  the  beam.  Both  contributions  to  the  total  frequency  spread, 

and  I  Aft  I  g ,  are  proportional  to  the  various  energy  spreads.  The  two 

_2 

different  proportionality  factors,  lal  in  the  case  of  lifting  and  y  in  the 
case  of  lAftL,  arise  because  the  intrinsic  particle  energy  spread  produces  a 

D 

rotational  frequency  spread  primarily  by  changing  the  particle's  radial 
position  whereas  the  various  energy  spreads  contributing  to  lAftig  merely 
result  in  a  longitudinal  velocity  spread.  Hence  the  various  contributions  to 
the  longitudinal  energy  spread  contribute  differently  to  the  frequency  spread. 

The  desired  stability  criterion  is  obtained  by  substituting  (6)  into  (5) 
and  becomes 

r  <  l'jic  [j  I  al  ( AE/E)  +  |(en/rt)2/2  -  C  /  rQ )  2y2  +  v/2y|/y2].  (9) 

Given  the  intrinsic  energy  spread,  beam  radius  and  emittance,  the  criterion  in 
(9)  implies  a  limiting  beam  current  which  if  exceeded  will  result  in 
instability.  As  a  simple  illustration  we  will  first  consider  the  negative 
mass  mode  in  a  low  current  conventional  betatron,  i.e.,  Ba  ■  0,  n  <<  1/2. 

o  S 
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From  (1),  the  dispersion  relation  for  the  negative  mass  instability 
2  2  2  2 

is  Aw.  ■  -  m  /oi ,  *  4i£  (v/y)(Z  (oi)/iZ  ),  where  we  have  assumed 

£  9  1  Ho 

that  I  Aw  I  «  la),  I  *  1//2.  Approximating  Z  (oj )  by  Z  (£w  )  we  find  that  the 

x  n  n  c 

growth  rate  is  r  ■  Z^Cv/y)1^2  I  (Z  (£o»c))  I ^/(iZ^^2.  Using  (9)  and 
neglecting  the  betatron  oscillations  we  recover  the  well  known  negative  mass 
stability  condition 


v  .  o  / AEn 2 

Y  -  IZ|(4«c)|  1  E ' 


(10) 


In  obtaining  (10)  we  have  assumed  that  the  real  and  imaginary  parts  of  Z^  are 
approximately  equal. 

Next  we  consider  the  full  dispersion  relation  (1)  for  arbitrary  b.  Here 

and  below  we  will  continue  to  neglect  the  effect  of  betatron  oscillations  on 

2 

the  stability  of  the  beam.  The  effect  is  generally  small  and,  due  to  the  b 
dependence  in  (6b),  it  favors  the  modified  betatron.  Therefore,  neglect  of 
the  betatron  oscillations  is  conservative  when  comparing  the  modified  and 
conventional  betatrons. 

If  we  may  continue  to  approximate  Z  (w)  by  Z  (£w  )  and  if  we 

U  j  r  '!  >  r  c 

set  5*0  then  (1)  is  in  general  a  sixth  order  (fourth  order,  if  b  ■  0) 
polynomial.  By  neglecting  resistive  wall  effects  we  are  assuming  that  these 
are  negligible  for  the  negative  mass  branch  of  the  dispersion  relation  in 
which  we  will  be  interested  here.  Of  the  six  roots,  two  pairs  correspond  to 
transverse  modes  (a  pair  being  a  forward  and  backward  wave,  essentially)  and 
one  pair  to  the  longitudinal  mode. 

For  a  given  set  of  parameters,  we  have  found  the  roots  numerically,  for  a 


large  range  of  £;  a  maximum  value  of  r/i  is  then  found  and  substituted  into 
the  stability  criterion  from  which  a  value  of  AE/E  is  computed.  An  important 


point  to  note  is  that  the  toroidal  field  has  more  of  a  stabilizing  effect  on 
the  high  l  modes  than  on  the  low  %  modes;  consequently,  as  b  is  increased  the 
behavior  of  the  impedance  for  small  l  tends  to  determine  the  maximum  V/l  and 
therefore  the  limiting  current. 

Typical  results  are  illustrated  in  Figs.  2  and  3.  Here  we  plot  the  beam 
current  vs  the  Lorentzian  full  width  energy  spread  required  for  stable  motion 
for  b  ■  0,  5,  10,  20,  rQ/a  »  6,  a/r^  -  3,  RQ/Z0  =  4,  and  Q  *  10.  Figures  2 
and  3  are  plotted  for  y  *  3  and  6  respectively.  The  dashed  lines  are  plots  of 


4- 
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For  consistency  we  must  be  to  the  left  of  these  lines.  This  restriction  may 
be  understood  by  considering  the  displacement  of  a  particle  from  the  center  of 
the  beam: 


6r  _  <SE/E 


+  betatron  oscillations 


where  5E  is  the  difference  between  the  particular  particle  energy  and  the 
average  energy  of  the  particles  in  the  beam.  It  follows  that 


fb  .  | AE/E  ■ 

r  -  l-2n  1 
o  s 


where  AE  is  the  full  width  of  the  distribution. 

The  effects  of  the  self  fields,  as  represented  by  ng  in  (6a)  and  (11)  are 
immediately  evident  in  the  plots  of  Figs.  2  and  3.  If  ng  were  zero  then  the 
dashed  boundary  lines  would  be  a  single  vertical  line  and  the  curves  would  be 


monotonic.  As  it  is,  the  effect  of  the  self  fields  is  basically  traceable  to 


the  increasing  (as  ng  increases  toward  1/2)  then  decreasing  (as  ng  increases 
beyond  1/2)  factor  lal*  The  multi-valuedness  of  the  curves  may  then  be 
understood  as  follows  (Refer  to  Fig.  4.):  For  very  small  currents  (Branch  I, 
in  Fig.  4)  the  cold  beam  growth  rate  is  small  and  an  increasing  function  of 
current.  The  self  field  index,  ng,  is  negligible  compared  to  1/2  and  the 
energy  spread  required  for  stability  is  an  increasing  function  of  current. 
This  is  the  regime  in  which  virtually  all  conventional  accelerators  operate. 
As  the  current  is  further  increased,  however,  a  second  branch  becomes 
accessible,  shown  as  3ranch  II,  in  Fig.  4:  While  the  cold  beam  growth  rate 
continues  to  increase  with  increasing  current,  the  intrinsic  energy  spread 
stabilization  mechanism  becomes  more  effective  due  to  the  increase 
in  lal  until,  near  ng  <  1/2  a  very  small  energy  spread  results  in  a  large 
spread  in  angular  velocity;  the  instability  is  therefore  easily  stabilized. 
More  simply  said,  we  have  stability  on  the  low  current  branch  (3ranch  I) 
because  the  growth  rate  is  small  and  on  the  high  current  branch  (Branch  II) 
because  the  stability  mechanism  is  strong. 

There  is  a  third  branch  which  appears  in  the  example  of  Fig.  2  above 
ng  *  1/2  (I  >  316  A)  for  b  ■  10  and  20  and  is  illustrated  as  Branch  III  in 
Fig.  4.  This  region  is  accessible  only  in  the  modified  betatron  since  we  are 
constrained  in  the  conventional  betatron  by  the  equilibrium  condition 
ng  <  1/2.  For  b  ■  5  in  the  example  shown  in  Fig.  2  the  stable  points  fall  to 
the  right  of  the  dashed  line  and  so  are  not  shown.  As  the  current  is 
increased  beyond  316  A  the  growth  rate  increases  and  the  stability  mechanism 
becomes  less  effective  as  lal  decreases.  Consequently  as  the  current  is 
increased  beyond  this  point  the  required  energy  spread  increases 
monotonically.  This  is  illustrated  by  the  third  branch  in  the  figure.  (The 
sharp  corner  between  branches  II  and  III  of  Fig.  4  will  probably  be  rounded 
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off  by  non-linear  terms  in  the  equations  of  motion;  such  terms  will  become 
important  for  n,.  *  1/2  since  this  condition  corresponds  to  the  vanishing  of 
linear  restoring  forces  on  the  beam  particles.) 

Finally  we  comment  on  the  accessibility  of  the  various  branches  available 
for  small  energy  spreads.  If  beam  injection  proceeds  slowly,  over  many  growth 
times,  say,  it  appears  that  only  the  lowest  branch  is  accessible;  attempting 
to  add  more  current  will  drive  the  system  unstable.  If  however  current  can  be 
introduced  into  the  accelerator  more  rapidly,  the  higher  current  branches  may 
become  accessible.  Only  a  carefully  designed  experiment  can  test  this 
speculation.  For  b  *  5  and  10  for  the  parameters  of  Fig.  2  typical  growth 
times  are  »  3  particle  circulation  periods  so  that  high  current  injection  on 
this  time  scale  is  a  practical  experimental  possibility.  The  third  branch  in 
Fig.  2  is  clearly  the  most  promising  for  very  high  current  operation. 

V.  Conclusions 

We  have  shown  that  the  addition  of  a  toroidal  magnetic  field  to  a 
conventional  betatron  may  significantly  improve  the  current  carrying  capacity 
of  the  betatron  by  controlling  the  collective  instabilities  which  limit  the 
current.  The  calculation  has  included  self  field  effects  and  a  simple,  though 
realistic  model  for  the  longitudinal  and  transverse  Impedances.  The 
stabilizing  effects  of  betatron  oscillations,  which  have  been  shown  to  include 
the  effects  of  emittance,  toroidal  geometry,  and  energy  shear  due  to  the 
electrostatic  potential  drop  across  the  beam,  become  stronger  as  the  toroidal 
field  is  increased,  given  a  fixed  beam  radius.  Inclusion  of  self  field 
effects  in  the  stability  criterion  (6)  has  been  shown  to  lead  to  a  multi¬ 
valuedness  in  tne  current  vs  energy  spread  plot  which  has  been  interpreted  as 
the  result  of  the  competition  between  the  growth  and  stabilization 
mechanisms.  Accessibility  of  the  high  current  branches  may  depend  on  the 
duration  of  the  Injection  process. 
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Figure  2.  Instability  threshold  current  vs  energy  spread  for  b  ■  0, 
10.  The  beam  and  chamber  parameters  are  rQ/a  ■  6,  a/r^  ■ 
3,  y  -  3,  Rq/Z  -  4,  Q  -  10. 
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Same  as  Fig.  2  except  y  m  6. 
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Figure  4.  Sample  plot  of  limiting  current  vs  energy  spread,  illustrating  the 


three  possible  branches 
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